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ABSTRACT 

Deep sub-millimetre observations using SCUBA are presented of the central regions of four 
high redshift clusters which have been extensively studied optically: CL0023+0423 (z = 
0.84), J0848+4453 (z = 1.27), CL1604+4304 (z = 0.90) and CL1604+4321 (z = 0.92). 
10 sub-millimetre sources are securely detected towards these four clusters at 850/im, with 
two further tenuous detections; the raw 850/im source counts exceed those determined from 
blank-field surveys by a factor ~ 3 — 4. In particular, towards CL1604+4304, 6 sources are 
detected with SgsouJy > 4mJy making this the richest sub-mm field discovered to date. Cor- 
rections for gravitational lensing by these high redshift clusters reduce these excess sources 
counts, but are unlikely to account for more than about half of the excess, with the remainder 
presumably directly associated with cluster galaxies. The 450 to 850/im flux density ratios of 
the detected sources are systematically higher (at a significance level > 98%) than those de- 
termined for blank-field selected sources, consistent with them being at the cluster redshifts. 
If subsequent identifications confirm cluster membership, these results will demonstrate that 
the optical Butcher-Oemler effect is also observed at sub-mm wavelengths. 

Key words: submillimetre — galaxies: clusters: general — galaxies: starburst 
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1 INTRODUCTION 

Our understanding of star formation at high redshift has under- 
gone a revolution in recent years. Lyman-dropout techniques (e.g. 
Steidel et al. 1996) have identified thousands of high redshift field 
galaxies with star formation rates of a few Mq yi~ , providing es- 
timates of the cosmic evolution of the mean global star formation 
rate (e.g. Madau et al. 1998). However, it has also become clear that 
dust plays a very important role in star-forming galaxies, absorbing 
a significant proportion of the emitted light. 

In order to directly address the role of dust, several sub-mm 
blank field surveys have been carried out using the Submillimetre 
Common-User Bolometer Array (SCUBA; Holland et al. 1999) on 
the James Clerk Maxwell Telescope (JCMT), to various depths and 
sky areas. These have provided n umber counts of su b-mm sources 
dow n to below 1 mJy at 850/im (jSmail et al. 1997|; Hughes et al. 
1998; |Barger et al. 1999j; plain et al. 1999j; |Eales et al. 2000t S cott 
et al. 2002; |Cowie et al. 2002| ), and have resolved the majority of 
the sub-mm backgro und radiation into ultraluminous galaxies with 
SFR > 100M©/yr ( frlain et al. 199S| ). The host galaxies of these 
sources show high optical obscuration, and so are missed by ultravi- 
olet selection techniques (cf. Adelberger & Steidel 2000). SCUBA 
observations pointed at Lyman-break galaxies have failed to de- 
tect the vast majority to a combined sample noise level of about 
0.5 mJy, indicating that these galaxies don't contribute significantly 
to the sub-mm background (e.g. Chapman et al. 2000). Using a sta- 
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tistical analysis, Peacock et al. ( 2000| ) demonstrated that Lyman- 
break galaxies are correlated with faint background fluctuations (at 
the ~ O.lmJy level) in a deep SCUBA image of the Hubble Deep 
Field, and calculated that on average the true star formation rates of 
Lyman-break galaxies are about a factor of 6 higher than those de- 
termined from their measured ultraviolet luminosities. This is con- 
sistent with other determinations that the global star formation rate 
at high redshifts should be at least a factor of a few higher than that 
measured optically (e.g. Hughes et al 1998). 

The ultraviolet and sub-mm selection techniques clearly sam- 
ple different populations of star-forming galaxies, and so to obtain 
a complete census of star formation at high redshifts, observations 
of both normal and dusty galaxies are required. Whilst these have 
been carried out for the field, to date high redshift cluster environ- 
ments have been neglected at sub-mm wavelengths. 

Rich galaxy clusters contain large numbers of galaxies at the 
same distance, and so provide important testbeds for models of 
galaxy evolution. The cores of low redshift clusters are dominated 
by a population of luminous early- type galaxies; the optical proper- 
ties of these evolve only very slowly with redshift, in a manner con- 
sistent with passive evolution of stellar populations which formed 
at a very early cosmic epoch (e.g. van Dokkum et al. 1998a,b, 
Schade, Barrientos & Lopez-Cruz 1997, Ellis et al. 1997, Stanford, 
Eisenhardt & Dickinson 1998). Low redshift clusters have been tar- 
getted at sub-mm wavelengths, using them as lenses to magnify the 
background population; only 2 cluster galaxies have been detected 
within 7 clusters at redshifts 0.19 < z < 0.41 (Edge et al. 1999; 
see also Chapman et al. 2002), and both of these were centre clus- 
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ter galaxies hosting powerful radio sources, in clusters with strong 
cooling flows. 

The composition of clus ters i s a strong function of redshift, 
however: Butcher & Oemler (1978) showed that at z > 0.3 a pop- 



ulation of bluer galaxies appears in many, but not all, clusters (the 
Butcher-Oemler effect). These blue galaxies often exhibit strong 
Ha or [Oil] 3727 emission lines, indicative of on-going star for- 
mation (e.g. Dressier et al. 1997, 1999). High redshift clusters also 
show an increased number of weak radio sources (e.g. Dwarakanath 
& Owen 1999, Morrison 1999, Smail et al. 1999, Best et al. 2002), 
whose presence is correlated with the Butcher-Oemler blue galaxy 
population, though without correspondence of the host galaxies. 
The host galaxies of these weak radio sources comprise a mixture 
of starbursts and AGN. This increase in activity in high redshift 
clusters is not unexpected. Low redshift clusters have the majority 
of their gas and galaxies in stable virialised orbits with the galaxies 
in the cluster centre having been stripped of gas (e.g. Gunn & Gott 
1972), switching off star formation. High redshift clusters, on the 
other hand, are much younger, still in their formation process with 
relatively high galaxy merger rates (e.g. van Dokkum et al. 1999) 
and a plentiful supply of disturbed gas, and thus provide ideal lab- 
oratories to induce both starbursts and AGN. 

To fully understand galaxy and cluster evolution, it is impor- 
tant to investigate the on-going star formation in these high— z 
clusters using sub-mm studies as well as the standard optical selec- 
tion techniques. This paper presents the first results of a sub-mm 
survey of the central regions of high-redshift clusters, designed to 
address this issue. The sample selection, observations and data re- 
duction are described in Section ^. In Section ^ the four clusters are 
individually discussed. The nature of the detected sub-mm sources, 
their likelihood of being cluster members, and their global proper- 
ties are discussed in Section^. Conclusions are drawn in Section^. 
Values for the cosmological parameters of Qm = 0.3, Qa = 0.7 
and Ho = 65kms _1 Mpc _1 are adopted throughout. 



2 OBSERVATIONS AND DATA REDUCTION 
2.1 Details of the observations 

Cluster targets were selected from the literature based upon the 
following criteria: (i) they should have redshifts z > 0.8; (ii) 
deep imaging of the field should have been carried out, confirm- 
ing a large overdensity of red galaxies; (iii) there should exist 
detailed spectroscopic observations of the field, with at least 10 
galaxies being spectroscopically confirmed as cluster members. 
These criteria should ensure that the cluster targets are indeed 
genuine high redshift clusters. Four clusters were selected accord- 
ing to these criteria and observing time constraints: CL0023+0423 
(z = 0.84), J0848+4453 (z = 1.27), CL1604+4304 (z = 0.90) 
and CL1604+4321 (z = 0.92). This sample of clusters is not in- 
tended to be complete, but should be representative. 

The central regions of the four cluster s were observed at 850 
and 450/im simultaneously using SCUBA ( Holland et al. 1999| ) on 
the JCMT during a variety of nights between February 2000 and 
December 2001 . Details of the observations are provided in Table [jj 
The 850/im observations cover an area of sky approximately 2.5 
arcminutes in diameter, and have an angular resolution of about 14 
arcsec. At 450 fim the observed sky area is slightly smaller and the 
beamwidth is about 7.5 arcsec. 

Observations were taken using the standard "jiggle-mode", 
whereby the secondary mirror is moved around 64-point hexag- 
onal pattern, with Is at each position, in order to fully sample both 



the 450 and 850/im image planes. The secondary mirror is chopped 
at 7 Hz between 'on-source' and 'off-source' locations, and the 
telescope is also nodded every 16 seconds so that the chop posi- 
tion is placed at the opposite side of the field; the field is therefore 
observed in an on-off-off-on pattern. A single 64-point jiggle se- 
quence takes 128 seconds (plus overheads) to complete. 

For this program the chop throw was fixed to be 45 arcsec in 
an east-west direction. This short chop throw results in the 'off 
beam positions often also being within the field of view; where this 
is the case, the sources will have associated negative counterparts of 
half intensity 45 arcsec east and west of the true source. In addition 
to indications that a short chop throw improves sky-subtraction, 
these negative regions provide an independent measurement of the 
source, increasing the signal-to-noise ratios and source reliabili- 
ties. The expected source density is sufficiently low that these re- 
gions are unlikely to coincide with another source. 

A blazar was observed approximately every hour through- 
out the observations, to correct the pointing of the telescope. The 
pointing measurements before and after each source observation 
were used to correct for any pointing shift during the observation. 
The sky opacities were monitored by regular skydip observations 
throughout the nights. Flux calibration was determined by observ- 
ing at least one of the flux calibrators Mars, Uranus and CRL618 
each night, using the same 45 arcsec chop-throw as the targets. It 
is estimated that the photometry at 850/im is accurate to about 5- 
10%; at 450/im the photometric accuracy is worse (predominantly 
systematics), with variations up to about 25%. 

2.2 Data reduction and cleaning 

The observations were reduc ed using the SCUBA software, SURF 
(Jenness & Lightfoot 199E). The reference measurements were 



subtracted from the signal beams and the bolometers were flat- 
fielded using the standard SCUBA flatfield. Where CSO Tau data 
was available and reliable (ie. stable, and consistent with the sky- 
dips) throughout the night, the 850 and 450/im opacities were 
calculated for each observa tion u sing the CSO Tau relations pro- 
vided by Archibald et al. ( |200C| ): r g5 o = 4.02(r CS o - 0.001) 
and T450 = 26.2(rcso — 0.014). Where this was not possible, 
opacities were calculated by interpolating between the values de- 
termined from the preceding and succeeding skydips. After extinc- 
tion corrections were applied, noisy bolometers, noisy integrations 
and strong spikes in individual bolometers were rejected; this re- 
moved between 5 and 10% of the data. The sky variations were 
removed by subtracting, for each bolometer at each time-step, the 
mean signal over all of the other bolometers excluding those which 
had significantly higher than average noise. 

The calibrated datasets from the different observing nights 
were weighted according to their rms noise levels and combined 
to produce maps of each cluster. These maps contain negative side- 
lobes east and west of each detected source, as discussed in Sec- 
tion [2T]. To remove these and to improve the fidelity of the images, 
the ma ps were cleaned using a modified version of the CLEAN algo- 
rithm (Hogbom 1974) commonly used in radio astronomy. In this 



method, first a normalised beam-map is created from observations 
of CRL618 with the same chop-throw. Then, the highest peak on 
the source map is found, and the beam-map scaled and subtracted 
at this position to remove 10% of this peak flux (ie. a loop gain 
of 0.1). This removes some of the positive detection of a source, 
and partially fills in its negative sidelobes. The highest peak is al- 
lowed to be either positive or negative: although negative values 
are unphysical, this is important otherwise overcleaning of the map 



© 2002 RAS, MNRAS 000, [[[{k] 



Sub-mm sources towards z ~ 1 clusters 3 



Table 1. Details of the SCUBA observations. For each cluster the redshift and pointing centre are given, together with the observation dates. For each 
observation date, the total number of 64-pt jiggle pattern repeats (10 jiggle repeats corresponds to about 30 minutes observing) and average CSO Tau (at 225 
GHz) are given. 



Cluster 



CL0023+0423 0.84 



RA Dec 
(J2000) 
00 23 52.3 +04 22 54 



J0848+4453 



1.27 08 48 34.6 +44 53 42 



CL1604+4304 0.90 16 04 25.1 +43 04 53 



CL1604+4321 0.92 16 04 31.5 +43 21 17 



Obs Date 
(dd-mm-yyyy) 

25- 10-2000 

27- 09-2001 
30-09-2001 
15-11-2001 
08-12-2001 

26- 02-2000 

27- 02-2000 

28- 02-2000 

02- 03-2000 

03- 03-2000 

04- 03-2000 

21- 03-2000 

22- 03-2000 
25-10-2000 

28- 01-2001 

29- 01-2001 

04- 05-2001 

05- 05-2001 

06- 05-2001 
08-05-2001 
27-09-2001 



N int 
(64-pt jigs 
130 
77 
50 
30 
90 
90 
78 
70 
65 
90 
96 
70 
70 
20 
70 
20 
40 
50 
35 
30 
61 



T CSO 

0.09 
0.06 
0.10 
0.06 
0.06 
0.07 
0.07 
0.05 
0.05 
0.04 
0.06 
0.04 
0.03 
0.09 
0.05 
0.05 
0.08 
0.07 
0.08 
0.08 
0.06 



can give rise to spurious features when only positive noise peaks 
and not negative ones are cleaned. This cleaning was repeated for 
a number (typically 50 to 100) of iterations, until the highest peak 
found was within 2.5a of the sky rms noise. 

Each residual map comprises sky noise plus sources fainter 
than the minimum cleaned flux density. The flux removed from 
this map consists of a set of delta functions around the original 
source peaks; these are then convolved with a Gaussian of the same 
resolution as the JCMT beam and added back on to the map. The 
combined maps were then smoothed slightly by convolving them 
with a Gaussian, to lower the resolution to 16.5 arcsec at 850/im 
and 9.5 arcsec at 450/im in order to enhance the detectability of 
faint sources. The 850^im and 450/im maps of the clusters so pro- 
duced are displayed in Figures |j] to Q Comparison with the initial 
raw maps confirms that the cleaning process has reduced the back- 
ground rms of the maps by removing the negative sidelobes, but not 
given rise to any new sources. 

2.3 Source extraction 

Perhaps the most uncertain aspect of the analysis of sub-millimetre 
data is the determination of the believability of sources from the 
final map. This arises primarily because of the small number of 
bolometers across the field (37 for the long wavelength array): 
noisy (or flagged) bolometers therefore lead to regions of higher 
noise in the map and can give rise to spurious sources. Several 
steps have been taken to avoid this. Firstly, all fields have been 
observed either close to transit or with some observations either 
side of transit so that the bolometers rotate around the field, thus 
minimising the effects of any bad bolometers. Secondly, except for 
J0848+4453, all fields have been observed during at least two sepa- 
rate runs, with consequent different noisy bolometers and different 
sky coverage. Thirdly, these different observing periods have been 
used to test the repeatability of detection of plausible sources (see 
below). Fourthly, on two of the 850^m maps (J0848+4453 at lo- 
cation 08 48 35.55, +44 54 24.7, and CL1604+4321 at location 16 



04 27.63, +43 21 52.9) the removal of bad bolometer data has led 
to regions of significantly shorter exposure and hence higher noise 
which may give rise to spurious sources: a 25 arcsec diameter re- 
gion has been masked around these two locations. 

From the final cleaned maps, both the peak and the integrated 
flux densities of likely sources were determined, the latter being 
measured by summing the flux density inside an aperture of ra- 
dius 15 (10) arcsec at 850 (450) /im. The peak flux density (in 
mjybeam -1 ) should be equal to the integrated flux density of the 
source (in mJy) for an unresolved source, and for extended sources 
it provides a lower limit. In low signal-to-noise observations, how- 
ever, the peak flux density can overestimate the true flux densities, 
and the integrated flux density usually provides a more accurate 
measure. Sources with both peak and integrated flux densities mea- 
sured at greater than 3.5a, where a is the rms sky noise of the map, 
were considered possible detections. Although within the 8 maps 
there are negative regions which reach —4a in peak flux density, 
in no case does a region reach —3.5a in both peak and integrated 
flux, and so there are unlikely to be any false positive detections at 
this significance level. 

For each field, the observations were then divided into two or 
three subsets; in the case of J0848+4453 this was done by separat- 
ing the first half of each night from the second half (ie. by bolome- 
ter location on the map), and in the other three cases the separation 
was made by observing date. Independent maps were made for each 
subset of data. These maps were then examined to investigate the 
reliability of the sources detected above. All of the sources detected 
above 4a in the overall map were also present above ~ 2a in each 
of the data subsets. 

The detected sources are listed in Table |[ The reliability of 
each source is also indicated in the table: those sources which are 
detected above 4a in integrated flux density are considered secure, 
whilst sources between 3.5 and 4a are less so. 
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Figure 1. SCUBA maps at 850 and 450/^m of the central regions of the CL0023+0423 cluster at redshift z = 0.84. Contour levels are plotted at (-4, -3, 3, 
4, 5, 6, 7, 8) times the sky rms noise level of each map: the 850 and 450/^m maps have rms noise levels of 1.2 mJy beam" 1 and 9.9 mJy beam -1 respectively 
(note that the maps have been smoothed). 



J0848+445 (z=1.27); 850 microns 
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Figure 2. SCUBA maps at 850 and 450/jm of the central regions of the J0848+4453 cluster at redshift z = 1.27. Contour levels are plotted at (-4, -3, 3, 4, 
5, 6, 7, 8) times the sky rms noise level of each map: the 850 and 450/im maps have rms noise levels of 1 .6 mJy beam - 1 and 12.2 mJy beam -1 respectively 
(note that the maps have been smoothed). A 25 arcsec diameter region centred on 08 48 35.55, +44 54 24.7 has been masked in the 850/tm map due to a higher 
noise level associated with bolometer flagging. 



3 DISCUSSION OF INDIVIDUAL SOURCES 



CL0023+0423 

This cluster was discovered opticall y and appears in the cluster 
sample of Gunn, Hoessel and Oke (19861. It forms part of the 



sampl e of 9 high redshift clusters which Oke, Postman and Lubin 
( 1998 ) have studied in detail with deep multi-waveband i maging 
and spectroscopic observations. Postman, Lubin and Oke ( 1998 ) 
show that the optical overdensity is actually the superposition of 



a cluster at z = 0.845, with a velocity dispersion of 415 km s 1 
(corresponding to a dynamical mass M ~ 3 x 10 14 M Q ) and 17 
spectroscopically confirmed galaxy members to date, and a group 
of lower velocity dispersion (a ~ 158 km s -1 ) at z = 0.827, with 
7 spectroscopically confirmed galaxies. These two structures are 
separated by less than 3000 km s -1 in redshift space, and may be 
interacting: about 57% of the cluster members show evidence for 



star formation activity (Postman et al. 199S ). 

The 850/^m SCUBA image of this cluster (Figure [l]) shows 
one source of 7.4 mJy integrated flux density (HzC10023_850.1), 
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Figure 3. SCUBA maps at 850 and 450/^m of the central regions of the CL1604+4304 cluster at redshift z = 0.90. Contour levels are plotted at (-4, -3, 3, 
4, 5, 6, 7, 8) times the sky rms noise level of each map: the 850 and 450^m maps have rms noise levels of 1.1 mJy beam -1 and 5.5 mJy beam -1 respectively 
(note that the maps have been smoothed). 
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Figure 4. SCUBA maps at 850 and 450/^m of the central regions of the CL1604+4321 cluster at redshift z = 0.92. Contour levels are plotted at (-4, -3, 3, 4, 
5, 6, 7, 8) times the sky rms noise level of each map: the 850 and 450/im maps have rms noise levels of 1.2 mJy beam -1 and 11.1 mJybeam - 1 respectively 
(note that the maps have been smoothed). A 25 arcsec diameter region centred on 16 04 27.63, +43 21 52.9 has been masked in the 850/tm map due to a higher 
noise level associated with bolometer flagging. 



which appears significantly extended; this source appears in all sub- 
sets of the data, although with positional offsets of up to 5 arcsec, 
presumably because of the extended nature of the source. A second 
source (HzC10023_850.2) is detected above 4er in peak flux, but 
below 4a in integrated flux, and so this is considered a marginal 
detection. No believable sources are detected at 450/xm. 



J0848+4453 

This cluster, the highest redshi ft in the current sample, was first 
discovered by Stanford et al. (1997) as a concentrated overden- 



sity of very red objects in a near-infrared field survey. There is a 
4.5a significance X-ray detection (L x ~ 10 44 ergs -1 , 0.5-2 keV) 
at the same location. Follow-up spectroscopy has, to date, con- 
firmed 10 galaxies as cluster members, with a velocity dispersion 
of 64 km s -1 suggest ing a cluster of about Abell richness class 



1 (Stanford et al. 1997). More recently, a second group of galax- 
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Table 2. The SCUBA sources detected from the maps, at each of 850 and 450/tm. For each cluster the map rms noise is provided (note that the maps have 
been smoothed), together with a listing of the RA, Dec, offset from the cluster centre, peak flux density, and integrated flux density of each source. Statistical 
eiTors on the flux densities of individual sources are the rms of the sky background for the peak flux densities and approximately 1 .5 times this value for the 
integrated flux densities. In addition there are associated calibration uncertainties for each map, at a level of 5-10% at 850/tm and up to 25% at 450/im. The 
penultimate column of the table indicates for each source whether the detection is fully secure: sources which are detected above 4cr in integrated flux density 
are considered secure, whilst sources between 3.5 and 4cr are less so. The 3.9<r 850/tm source towards CL1604+4304 is also classified as secure because it is 
detected at 450/im. The final column gives the 450/850/im flux density ratio (or upper limit) for 850/tm sources with secure detections which lie within the 
450/im field of view. 
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i es has been detec ted only 4 arcminutes away, with redshift 1.26 
(Rosati et al. 19991. This group also has associated X-ray emission 



of similar luminosity, and spectroscopic redshifts have confirmed 
cluster membership for 6 galaxies to date. The J0848+4453 cluster 
therefore appears to be part of a massive high redshift structure. 

The SCUBA maps of the central regions of this cluster indi- 
cate thepresence of two luminous sub-mm sources at 850/im (see 
Figure g; Table ^). A slightly extended source is seen to the west 
of centre of the map (HzC10848_850.2), and a very luminous (> 
13mJy) source close to the north-eastern edge (HzC10848_850. 1). 
This latter source is amongst the most luminous 850/tm sources 
known; although it lies close to the edge of the field, it occurs in all 
subsets of the data and so its detection is considered secure. 



CL1604+4304 

Like CL0023+0423, the cluster CL1604+4304 (z = 0.9 0) was first 
discovered in the optical imaging search of Gunn et al (19861 and 



subsequently formed part o f the sample of 9 high-redshift clusters 
studied by Oke et al ( 199S ). This c luster was studied in consider- 



able d etail by Postman et al (1998) and Postman, Lubin and Oke 



(2001 ): spectroscopic redshifts have been obtained for over a hun- 



dred galaxies in the field, with 22 being confirmed cluster mem- 
bers. The velocity dispersion of these galaxies, a ~ 1220 km s _1 , 
is consistent with an Abell 2 richness cluster although the clus- 
ter is only marginally detected in X-rays with a luminosity of 
~ lO ^erg s" 1 (0.1 to 2.4 keV; Castander et al. 1994), and Smail 



et al. (1994) failed to detect a weak lensing signal towards it. Ap- 



proximately half of the confirmed cluster galaxies are emission line 
(star forming) objects. 

This cluster is particularly interesting as it forms part of one 
of the most massive structures known at high redshifts. The equally 
rich cluster CL16044-4321 (also in the SCUBA sample: see below) 



is offset only 17 arcmins on the sky and 4300 km s -1 in redshift, 
and there is a large overdensity of red galaxies in a 'filament' join- 
ing the tw o clusters, includin g another smaller condensation of red 
galaxies (Lubin et al. 2000). Together these seem to comprise a 



very massive high-redshift supercluster. 

The 850/im SCUBA map of the central regions of this cluster 
is quite remarkable (Figure ^). This is both the deepest map in the 
sample, and the one taken in the best, most stable, sky conditions 
and therefore of high image integrity. Five sources are clearly and 
securely detected with flux densities in excess of 4.4 mjy (4cr; see 
Table |). A sixth source (HzC11604a_850.6) is detected slightly be- 
low 4<r significance, but appears as a luminous source in the 450/im 
observations (HzC11604a_450.3), confirming its reality. A further 
two of the 850/im sources are also detected at 450/im, providing 
still further confirmation of their reality. A fourth 450/im source 
is marginally detected at between 3.5 and 4a, but the disparity be- 
tween its peak and integrated flux densities and its non-detection at 
850/im cast doubt on its reality. The rms of the 450/iJy map is sig- 
nificantly lower than those of the other clusters in the sample due 
to the low opacities during these observations. 

The detection of 6 sources brighter than 4mJy at 850/im 
makes this the richest sub-mm field observed to date. The source 
density on this map is a factor of 5 higher than expected from blank 
field counts (see Section 1.1). 



CL1604+4321 

As discussed above, CL 1604+4321 at redshift 0.92 forms part of 
a potential hig h redsh ift supercluster together with CL 1604+4304. 
Postman et al. ( 2001 ) studied CL 1604+4321 in detail, and spectro- 
scopically confirmed about 40 galaxies as being associated with the 
cluster, with a velocity dispersion of 935 km s . 

One source is clearly detected above the 4a level at both 850 
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and 450/im in the SCUBA maps (HzCl 1 604b_850. 1 and 450. HfTF 
ure ^; Table |^). A second source, towards the north-east edge of 
the 850/im map, is detected with an integrated flux densit} level 
between 3.5 and 4a, and so is considered to be an insecure detec- 
tion. 



4 RESULTS 



Ten sub-mm sources have been securely detected towards 
four high redshift clusters at 850/im, with a further two 
detections. The important issue is now to determine whether 
sources are indeed associated with the clusters or whethe 
are merely projected along the same line of sight. Before 
the sources are identified and their redshifts measured, this 
tion cannot be unambiguously answered, but in this section 
indications are given based upon the source counts and 
450 / 850/im flux density ratio which serves as a rough 
dicator. 
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4.1 Integrated Source counts 

The integrated source counts at 850/im and 450/im have 
termined as a function of limiting flux de nsity u sing the i 
flux densities of each source. Scott et al. (2002) found, usin: 



ulated maps, that the flux densities of sources in the 8 mJy 
are typically overestimated due the effects of both noise and 
confusion; here no corrections were made for this effect, 
our sky area covered is too small and contains too many soui 
accurately calculate any such corrections. At the significance 
(> 4cr) of the sources considered, the flux boosting effects 
noise should be minimal, and in any case any such flux bo 
is compensated by the fact that the 15 arcsec radius apertme 
not contain the total source flux density — in theory only 
for a point source. Confusion effects, however, could be 
if faint sources below the flux limit combine to give an 
source above that limit. One or two of the more extended 
could conceivably be the result of so urce confusion, althou; 
simulations of Hughes & Gaztanaga ( |200C| ) suggest that this 
is of little importance at flux densities above ~ 2.5 mJy, i 
from here on it is discounted. 

Only the securely detected sources were considered 
analysis. At 850/im the sky area observed in each map was 
to be 4.5 square arcminutes down to a flux density level of 4 
the rms noise quoted in Table ^, plus a further 1 square 
around the edges of the map to a flux density level 50% 
At 450/im the sky areas were similarly 3 square arcminutes 
highest sensitivity and 1 square arcminute at the lower sens: 
Between each of the 8 limiting flux densities (that is, the 
values for the central and edge regions of each map), the 
sampled brighter than this limit was determine and used to cal 
the number density of sources within that flux density range, 
these were calculated the number density of sources above 
flux density levels at 850/im and a single level at 450/im 
provided in Table 3. It is important to note that at this stage 
count is taken of any possible gravitational lensing by the 
and so the raw source counts likely overestimate the true c 



sk/ 



this is discussed in detail in Section 4.2 . 

The integrated source counts determined at 850/im arel com- 
pared with those determined from blank field surveys on Fi£ ure |^. 
The comparison data are taken from the Hubble Deep Field ; urvey 
( Hughes et al. 199E ), the various SCUBA lensing surveys (|Smail 
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Table 3. The integrated source counts as a function of limiting flux density 
at 850 arjd 450/im as observed in the fields of the high redshift clusters. 
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Blain et al. 19% ; "owie et al. 2002), the Hawaii Sur 



are clearly in excess, by a factor -~3-4, of those in the 
surveys. Note that a similar excess of 850/im counts has 
found in the field surrounding the high redshift (z — 3.8) 
laxy 4C14. 17, which is thought to reflect a high-redshift 
uster ( Ivison et al. 2000 ). It should be noted that the sky ar- 



are rela tively small, and SCU BA sources are known to 
clustered (Almaini et al. 2002), which could give rise to 



lan average counts being measured in some locations. Al- 
I his possibility cannot be definitively excluded until identifi- 
md redshifts are available for the sources, for such a strong 
ity to occur by chance just at the location of these clusters 
e an odd coincidence, and so it is far more likely that the 
real (or due to lensing; see below). 

450/im determining the source counts is considerably more 
due to the fewer sources d etecte d and the larger uncer- 
in calibration. Smail et al. (2002) derive 450/im source 
be 2100 ± 1200 deg -2 sources brighter than lOmJy and 
00 deg" 2 sources brighter than 25 mJy. The 3800 ± 1900 
2 determined here to the same 25 mJy flux density limit 
an order of magnitude in excess of their counts, although 
s are clearly large. 



ceg 



4.2 G - avitational lensing? 



these sources are observed in directions towards concen- 
of mass, the possibility of gravitational lensing by the high 
clusters must be considered. Gravitational lensing has two 
effects upon the observed distribution of source counts: 
luminosity of each lensed source is increased, boosting 
rent source counts, but the sky area observed is decreased 
regions. 

to any lensing the total source counts observed will be 
J* max n(S u , z)dz, where n{S v , z) describes the source 
a function of redshift. After lensing by a factor /i(#, z), 
is the angular offset from the lens centre, the observed 
ounts will be 



le used 
ipr 

s 



the 



S v 
be 



dz. 



source counts can be described by a power law, 
S„ a , then n( — , z) = n a n(S v ,z) and therefore 
9) = ii(6) a - 1 N(S^). When the slope (a) of the power- 
unity then the effects of gravitational lensing can boost 
source counts. Figure ^ demonstrates that the sub-mm 
ounts indeed have a very steep slope: the slope of the inte- 
spurce counts, N(> S v ), is about 1.8, and so a ~ 2.8. Thus 
of a factor 2 provides a boost in the source counts 
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Figure 5. The integrated source counts at as observed in the fields of the high redshift clusters (filled circles) compared with those determined from 

various SCUBA field surveys (open symbols), taken from the references indicated in the lower left corner of the plot. A clear excess of sub-mm galaxies is 
seen in the direction of the high redshift clusters. 



of approximately a factor 3.5. It is this dramatic effect which has 
led to low redshift rich clusters being used as lenses to study the 
faint sub-mm source population (e.g. Smail et al 1997). 

Modelling the lensing object as a singular isothermal sphere 
(SIS), that is, one with a central density profile of p(r) oc r~ 2 , the 
factor by which the sources are magnified at an angular offset 8 is 
related to the Einstein radius, 8e, by: 



i 



The Einstein radius can be expressed in terms of the properties of 
the source and lens as (e.g. Blandford and Narayan 1992): 

47TCT 2 £> LS 

c 2 Dos 

where a is the velocity dispersion of the lens, Dls is the angu- 
lar distance between the source and lens, and Dos is the angular 
distance between the source and the observer. The lensing mag- 
nification is maximal for sources observed at the Einstein radius 
(theoretically infinite for a point source, and in practice giving rise 
to strongly-lensed arcs or multiple images) but falls off rapidly be- 
yond this: at 8 = 2#e the magnification factor p(8) ~ 2, and by 
4#e it is only a factor of 1.33. For flatter central density profiles, 
for example the NFW profile (Navarro, Frenk & White 1997) with 



p(r) 



the magnification factors beyond 8e are lower than 



for the SIS. 

The velocity dispersions of the four clusters in this sample are 
very uncertain because of the small numbers of redshifts currently 
available. The est imated values rang e from 415 to 1220 km s' 1 for 
the four clusters (Rosati et al. 1995; Postman et al. 2001). Adopt- 
ing a velocity dispersion of 1000 km s - , a lens redshift zl = 1 
and a source redshift zs = 3, the Einstein radius corresponds to 



#e ~ 13arcsec. This value is significantly smaller than for low 
redshift rich clusters (up to 30 or 40 arcsec), firstly because the 
high redshift clusters are of significantly lower mass (velocity dis- 
persion) and secondly because the optimal lensing configuration 
has equal source-lens and lens-observer distances which favours 
0.2 < zl J; 0.4 for 2 < zs < 4. Therefore the effects of magnifi- 
cation by these high redshift clusters are of lesser importance. 

A proper determination of the unlensed source counts can- 
not be carried out without an accurate mass model for these clus- 
ters, but a rough judgement of the importance of lensing magnifi- 
cation can be made. Integrating the expression for the magnified 
source counts across the field using an Einstein radius of 13 arcsec 
(a ~ 1000 km s ) predicts a total boosting of the 850/im source 
counts of a factor between 1.5 and 2 for an SIS profile. This is 
consistent with the boosting in number counts that might be ex- 
pected by examining the distribution of detected sources within 
the fields: the offsets of each individual source from the cluster 
centre are given in Table ^. HzC11604a_850. 1 has a angular off- 
set from the cluster centre of only ~ 20 arcsec indicating that it 
is likely to be significantly magnified; the relatively low 450 to 
850/im flux density ratio of this source (see Section |4.3j ) further 
supports the suggestion that this is a background source. Four fur- 
ther objects (HzC10023.850.1, HzC10848.850.2, HzC11604a.850.3 
and HzC11604a_850.4) have angular offsets of between 30 and 40 
arcsec, with consequent flux boosting of ~ 50% if these are lensed 
background objects, but all of the remaining 7 sources are more 
than 40 arcsec away from their cluster centres, with expected flux 
boostings of < 25%. The radial distance of all but one sources is 
at least double the Einstein radius, making it extremely unlikely 
that any multiple image strong lensing affects the number counts in 
these clusters. 

It should be stressed that, unless the current velocity disper- 
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sion measurements greatly underestimate the mass of these clus- 
ters, these calculations represent the upper limit of the boosting 
that gravitational lensing by the cluster could feasibly provide. The 
a ~ 1000 km s adopted is very much at the maximum for clus- 
ters at these redshifts, and the SIS profile represents the 'worst- 
case' scenario. Further, for CL 1604+4304 with the most detected 
sub-mm sources and the highest estimated velocity dispersion, not 
even a weak-lensing si gnal w as seen for optical galaxies in the 
analysis of Smail et al. (1994). Therefore although cluster gravi- 



tational lensing is likely to have a significant effect upon the source 
counts, this is not nearly to the extent of the factor ~ 3.5 higher 
source counts that is observed. 

An alternative factor that could conceivable boost the source 
counts is gravitational lensing by individual galaxies. It has been 
suggested that the brightest SCUBA source in the Hubble Deep 
Fiel d is gravitationally lensed by a z ~ 1 elliptical (D ownes et al. 
1999), and within a z ~ 1 cluster there is a higher surface den- 
sity of elliptical galaxies which could act as lenses. The Einstein 
radius of an individual elliptical galaxy at these redshifts will be 
only #e ~ 1-2 arcsec; an examin ation of optical data available for 
CL1604+4304 (|Smail et al. 1994 suggests that the SCUBA source 



lies this close to a bright elliptical in only one of the six cases, 
that being HzC11604a_850. 1 which was already suggested to suffer 
strong cluster lensing. Once again, therefore, this effect is unlikely 
to significantly affect the number count statistics, but could be rele- 
vant in isolated individual cases. This is an important effect to bear 
in mind during follow-up studies as it can easily lead to misidenti- 
fication of the host galaxy (cf Downes et al.1999). 

All in all, the conclusion must be (particularly when coupled 
with the redshift indicators below) that the higher source counts are 
in part due to lensing effects, both from the cluster and individual 
galaxies, but that a significant population of cluster sources is also 
detected. 



4.3 Redshift indicators 

For a dust temperature of 30 to 50K, the grey-body dust spec- 
tral energy distribution (SED) peaks at ~ 100/itri. Thus, whilst 
the 850/im data point lies on the Rayleigh-Jeans tail, providing a 
strong positive K-correction which results in the 850/iin flux den- 
sity of a given source appearing the same at all redshifts 1 < z ;$ 
5 — 10 (e.g. Blain & Longair 1993), the 450/^m data point moves 
steadily through the peak of the SED with increasing redshift. The 
450 to 850/im flux density ratio therefore decreases monotonically 
with redshift and so can be used as a broad redshift indicator. 

A note of caution needs to be added when applying this re- 
lation in practice, since 450/im calibration is notoriously difficult 
except in the best atmospheric conditions (such as those for the ob- 
servations of CL1604+4304), and calibration errors would lead to 
systematic offsets, and hence uncertain results. The radio to sub- 
mm spectral index (e.g. Carilli & Yun 2000) is a better tool to use 
as a redshift estimator, and in addition the accurate positions deter- 
mined from the radio data allow optical identifications to be made. 
Radio observations of these fields are proceeding and will provide 
more definitive results in the future, but with the current data the 
450 to 850/im flux density ratio can provide some good first indi- 
cations. 



Lutz et al. (2001) investigated the how the 450 to 850^m flux 



density ratio would vary as a function of redshift from the spec- 
tral energy distributions of 14 local ultra-luminous infrared galax- 
ies (ULIRGs) with existing far-infrared and sub-mm photometry 
(from Klaas et al. 2001). The relation they derived is shown in 



Figure ^a. 450 to 850/^m flux density ratios were calculated for 
all of the securely detected 850/xm sources in the current sample 
within the sky area of the 450/xm maps. These ratios are tabulated 
in Table ^ and are displayed as a histogram in Figure ^b. For 4 
sources a flux density ratio is determined, for 3 of which the ratio 
is high, > 5, consistent with being at redshifts z < 1. A further 4 
sources have determined upper limits on their flux ratios, which in 
3 cases suggest relatively high redshifts when compared directly to 
the predictions of Figure ^a, although the presence of any cooler 
(T ~ 10 — 15 K) dust components within the galaxies would give 
rise to lower 450 to 850/im flux density ratios at lower redshifts. 

In Figure ^c a similar histogram of flux density ratios has been 
derived for sources detec t ed in blank field surveys, using the 8 mjy 
survey ( Scott et al. 2002 ; Fox et al. 2002 ) and the SCUBA lensing 
survey (e.g. Smail et al 1997, 2002); these provide a good spread of 
sources both more and less luminous than those of the high redshift 
cluster sample, and so provide a good comparison sample. This fig- 
ure provides a strong visual indication that the sources in the cur- 
rent sample have systematically higher 450 to 850/im flux density 
ratios than those in the blank field surveys, for example considering 
the fraction of sources with 450 to 850^im flux density ratios above 
and below 4.0 in each sample. Survival analysis (e.g. Feigelson & 
Nelson 1985) can be used to statistically investigate datasets such 
as these with large numbers of upper limits. The Kaplan-Meier es- 
timator gives a mean 450 to 850/im flux density ratio for the blank 
field selected sources of 3.0 ±0.2 and even if all of the 5 upper lim- 
its for the cluster sample are between 2.0 and 2.5Q then the mean 
for that s ample will be > 4. U sing the ASURV survival analysis 
package (La Valley et al. 1992) the likelihood of the sources from 



the cluster and blank field surveys being drawn from the same par- 
ent population, as determined from the 450 to 850/im flux density 
ratio distributions, was calculated according to 5 different survival 
analysis methods based around a generalised Wilcoxin test: taking 
the mean of these results, the probability of the parent samples be- 
ing the same is < 1.8%. 

To summarise, the 450 to 850pim flux density ratios of the 
current sample are systematically higher than those of blank field 
surveys. This is most likely due to an extra population of sources 
associated with the high redshift clusters, which have higher 450 to 
850/iin flux because their redshifts are significantly lower (z ~ 1) 
than sources typically detected in blank field surveys (z ~> 2; e.g. 
Smail et al. 2000, Fox et al. 2002). 



4.4 Dust masses and star-formation rates 

The mass of warm dust, Md, can be calculated from the sub- 
millimetre flux using the equation: 



M d 



s (u b S ) dI 



(1 + z)Kd(l , Test)B(u 

rest , J d I 



where S is the observed flux density, u b s and !/ rcst are the ob- 
served and rest-frame frequencies, Z?l is the luminosity distance, 
z is the redshift, Kd is the mass absorption coefficient, B is the 
black-body Planck function, and Td is the dust grain temperature. 
Here, a dust temperature Td ~ 40 K, and a mass absorption co- 
efficient K d = 0.067(i/ r cst/250GHz) /3 m 2 kg" 1 with f3 = 2 are 
adopted: for a discussion of the uncertainties on these parameters 
see Hughes, Rawlings and Dunlop (1997) . With these properties 

T The Kaplan-Meier estimator cannot be applied to the cluster sample since 
the lowest data point is an upper limit rather than a measured value. 
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then a source with a flux density of 10 mJy at 850/im would have a 
dust mass of ~ 5 x 10 s Mq if it were at redshift z ~ 1. If it is as- 
sumed that this dust is heated primarily by young stars, then scaling 
the sub-mm luminosity to the star-formati on rate by using nearby 
starbursts such as M82 (Hughes et al. 1997), the star formation rate 
for a lOmJy 850/im source at z = 1 is ~ 2OOOM0yr _1 . Thus 
any confirmed sub-mm emitting cluster members are undergoing 
extremely active starbursts. 



5 CONCLUSIONS 

SCUBA observations have been presented of the central regions of 
four rich cluster of galaxies at redshifts z ~ 1. Compared with the 
number counts determined from blank-field surveys, an excess of 
sub-mm sources of a factor ~ 3 — 4 at and up to an order 

of magnitude at 450/xm has been detected. Some fraction of this 
excess will be due to the effects of gravitational lensing magnifi- 
cation of background sources by the clusters; an analysis of pos- 
sible source count boosting by gravitational lensing has been car- 
ried out and demonstrates that, unless the measured velocity dis- 
persions greatly underestimate the cluster masses, this could boost 
the source counts by at most a factor of two. The residual excess 
counts are most likely directly associated with the high redshift 
clusters. This conclusion is strengthened by the typically higher 
450 to 850/iin flux density ratios of these sources as compared with 
blank-field selected sources. These higher ratios are consistent with 
a significant proportion of the sources being at lower redshifts than 
those of the blank field surveys (typically z > 2). 

Confirmation that these sub-mm sources are indeed asso- 
ciated with cluster galaxies would indicate a dramatic increase 
with redshift in the proportion of strongly starbursting galaxies 
within clusters. This mirrors the increasing fraction of blue, star- 
forming galaxies indicated by the Butcher-Oemler effect (Butcher 

6 Oemler 1978), and the recent discovery that more distant clus- 
t ers often contain an enhance d population of weak radio source s 
(bwarakanath & Owen 19991; Morrison 19991; best et al. 20021). 



All of these results demonstrate that high redshift clusters are still 
in a very active state of formation. With a larger dataset it would 
be particularly interesting to investigate the radial distribution of 
the sub-mm sources within the clusters, and to investigate whether 
mergers between clusters, which are common at these high red- 
shifts, increase the proportion of starbursting galaxies. 
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Figure 6. Top: The 450 to 850/tm flux density ratio as a function of redshift 
derived by Lutz et al. (2001) by redshifting the spectral energy distributions 
of 14 local ultra-luminous infrared galaxies (ULIRGs). The solid line dis- 
plays the median relation and the shaded region shows the total scatter of 
the points. Middle: Histograms of the 450 to 850^tm flux density ratios 
for those sources detected with a significance above 4cr in the fields of the 
high redshift clusters. The filled boxes represent sources detected at both 
wavelengths and the open boxes represent the 4cr upper limits for sources 
undetected at 450/^m. Only the measured ratios or limits are shown; statis- 
tical and calibration errors (which can be > 25%) are not indicated. Bot- 
tom: A similar plot combined for sources in the blank fields of the 8 mJy 
survey (Scott et al 2002; Fox et al 2002) and the SCUBA lensing survey 
(Smail et al 1997, 2002). The sources detected towards the high redshift 
clusters have systematically higher 450 to 850^tm flux density ratios, con- 
sistent with being at lower redshifts (z ~ 1) than those found in blank field 
surveys (typically z > 2). 
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